Dynamics of wave packets for the nonlinear Schrodinger equation in random potential 
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The dynamics of an initially localized Anderson mode is studied in the framework of the nonlinear 
Schrodinger equation in the presence of disorder. It is shown that the dynamics can be described in 
the framework of the Liouville operator. An analytical expression for a wave function of the initial 
time dynamics is found by a perturbation approach. As follows from a perturbative solution the 
initially localized wave function remains localized. At asymptotically large times the dynamics can 
be described qualitatively in the framework of a phenomenological probabilistic approach by means 
of a probability distribution function. It is shown that the probability distribution function may be 
governed by the fractional Fokker-Planck equation and corresponds to subdiffusion. 
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In this work the dynamics of an initially localized An- 
derson mode is considered. It is described by the nonlin- 
ear Schrodinger equation (NLSE) in the presence of dis- 
order [l|, [3, [J. In the linear one-dimensional case, eigen- 
functions are localized This problem is relevant to 

experiments in nonlinear optics, for example disordered 
photonic lattices @, 0] , where Anderson localization was 
found in the presence of nonlinear effects, as well as to 
experiments on Bose-Einstein Condensates in disordered 
optical lattices |8|, |9|, |lOi lUl, IjJ ■ It was shown that the 
presence of nonlinearity leads to an essential complication 
of a mechanism of localization [l|, [3, [3, > and the inter- 
play between disorder and nonlinear effects leads to new 
interesting physics [1, [M Q E H [13, E Hi!, HI . 
In particular, the problem of the spreading of wave pack- 
ets and transmission are not simply related [23L [24 | , in 
contrast with the linear case. In spite of the extensive 
research, many fundamental problems are still open. In 
particular, the spreading of an initially localized wave 
packet in nonlinear finite chains was challenged in numer- 
ical studies with realizations of subdiffusion [l^, [ItI, Ull 
and discrete breathers [H, [13, [2l|, [H, [2^ . The absence 
of the wave packet diffusion was observed as well and 
this behavior of initially localized wave packets was ex- 
plained by quasiperiodic solutions in the long time limit 
dynamics |25]. 

The system under consideration is the NLSE 



(1) 



where is a nonlinearity parameter. The variables are 
chosen in dimensionless units and the Planck constant is 
h = 1. The random potential V = V{x), x € (—00, -l-oo) 
is such that for the linear case (/? — 0) the Ander- 
son localization takes place, and the system is described 
by the exponentially localized Anderson modes (AM)s 
^t^ij^^ = vE'fc(a;), where ^'^^^(a;) are real functions and the 
eigenspectrum cjk is discrete and dense Therefore, 
the problem in question is an evolution of an initially lo- 
calized wave function il;{t = 0) = 'i>ig{x). Projecting Eq. 
([1]) on the basis of the Anderson modes 

^Pix, t)^Y. ^ E Ckmkix) , (2) 



we obtain a system of equations for coefficients of the 
expansion Ck 

idtCk = oJkCk + !3 Yl M^)ClCk,Ck, , (3) 

ki ,k2,k3 

where A(k) is an overlapping integral 

yl(k) = A{k,ki,k2,k3) 
*fc(a;)*fe,(a;)*fe,(a;)^'fe3(a;)da;. (4) 

The initial conditions for the system of Eqs. ([3]) are 

Cfc(t = 0) = a^^ = flfe = Sk,io ■ (5) 

Equations ([3]) correspond to a system of interacting non- 
linear oscillators with the Hamiltonian 



(6) 



Therefore, Eqs. ([3]) are produced by the Poisson brackets 
{if, . . .} pB by means of the Liouville operator 



f-lsH -If^ _d__dH_ _d_ 



(...). 

(7) 

Since LH = and i7({Cfe,C^}) = H{{ak,al}), we ob- 
tain that the Liouville operator is an operator function 
of the initial values: 



dHiak,al) d 9H(afc,a*) 



da.*. 



dak 



d 
d^k 



(8) 



and corresponds to the linear equation dtC = LC. Thus, 
the Liouville operator is the following combination: 



L ^ -i{Lo + (3Li) , 
where Lq = J2k (afegf^ ~ c.c.^ and 

d 



(9) 



dak^ 



c.c. 
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Here c.c. denotes a complex conjugation. Finally, we 
obtain that the system of nonlinear ordinary differential 
equations ([3]) is replaced by a system of linear partial 
differential equations: 

dtCk{t) = LCk{t) , fc = l,2,...,. (10) 

A formal solution of Eq. (fTi]) is the expansion 



ji=0 



ai,=Sk 



(11) 



The nonzero contribution to the first power over t of the 
expansion pT|) is due to the term 



' dak 



(12) 



while (Li — L^^^)ak = is due to the initial conditions 
([5]). Moreover, the contribution of Li — i^*'-' without L^^'* 
is zero in all powers of the expansion pip . For example, 
the nth power term for k ^ Iq is 



^ A{lo,lo,l,l)\ai^^\'^d^ 



ak = i''A^{lo,lo,k,k)Sk,io 



It has a non zero contribution only in the power of the n+ 
1 order after the action of the L^^^ term [23|. Therefore, 
keeping only the Lf'^ term in Eq. p4p means neglecting 
0{(3'^t^) terms in the expansion (jll[) . This solution is 
valid up to a time scale t < = 1/ [3. 

To obtain a solution in the framework of this approx- 
imation, first we eliminate the linear term Lq from Eq. 
(fTO| by substitution 



(13) 
(14) 



Cfc(0 =exp(-^Lo^)Cfc(^)• 
After this substitution, Eq. (jTU]) reads 

Taking into account that 

exp[— iLo^] — exp 



is the phase shift operator for the complex values au — 
|a/c|e*"^'=, we obtain 



Li(t) = ^A(k) 



exp[iAwt]o^^afc2afc3 
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(15) 

where Aw = t^^fe^ + lou^ - ujk^ - cj^a ■ 

Solutions of Eq. (fT4|) for k ^ Iq are functions which 
are zero at t = 0. These are 



Ck (t) ^ ak + 



-i0Ao\a,g\^t _ iAujt 



Here Aq = A{Io,Iq,Io,Io) and Ai = Ai{k) = 
A{Iq, lo,la,k) , k Iq, while Aw now is Aw = ujk — oji^. 
The complex conjugation of Eq. (fT6|l is a solution as well. 
A solution for k — Iq is a function of <pig — l3Ao\aig\'^t, 
which corresponds to the initial conditions Eq. ([5]): 



Cioit) = aig exp(-j/3Ao|a;„pt) . 



(17) 



Using these analytical solutions for the coefficients Ck (t) 
and Eq. (|13p . one obtains the solution of the of NLSE 
([T]) in the first order approximation over 1//3 as a sum 



tpit) = aig exp{-iuJi^t)'i'ig{x) - A(3\aig\^ai„ 

S 

^Ek'Mk)- 



*fe(a:), (18) 



where lui^ — ujig +f3Ao\a\'^ and wi^ = oJif, — f3AQ\a\'^, while 
prime means that k ^ Iq. When /3 — > 0, we have at the 
asymptotically large times tp ^ oo that uji-^ = uii^ — ojiq, 
and the sine function is 



2 sin 



lim 

t~*oo UJk — UJI2 



The sum in Eq. p8)l equals zero. Therefore, for /? = 0, 
one obtains '0(i) = e~^'^'^o'i>ig{x) that corresponds to a 
solution of the linear problem. 

For nonzero values /3 and t < tp the sine function can 
be approximated by tp for LOk ~ w;^. Then summation 
in Eq. (fT5)) can be estimated by adding and subtracting 
the term with k = Iq. Using the definition of the over- 
lapping integrals Ai{k) and J2k *fc(a;)*fc(2/) = S{y- x), 
one obtains an approximation for Eq. (|18p 

V^i) ^ vl/;^(x)e-*"'i* - 4/3t[vE'3 (x) - Ao^i„{x)] sin(w,,0 . 

(19) 

It means that at t < tp the wave function is localized 
and its evolution corresponds to the periodic oscillations 
with the frequencies and uJi^. It is worth mentioning 
that Eq. (|18p is valid for both finite and infinite systems 
(either discrete or continuous). 

Consideration of the dynamics beyond t > tp relates 
to the calculation of nonzero contributions of operators 



Li 



■-(0) 



and 



^(0) 



acting on the initial condi- 
tions. This combinatorics leads to essential difficulties 
for analytical treatment. To overcome this obstacle the 
dynamics of the initially localized states can be consid- 
ered qualitatively in the framework of a phenomenolog- 
ical probabilistic approach. To explain how the prob- 
abilistic approach works, let us demonstrate it first for 
the localized solution of Eq. . One obtains from the 
expansion (jlip 



00 
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n=l 
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(0) 



(20) 



The operator L^^^ corresponds to the population of all 
(16) states 4'^ by transitions from the state '^i^. Since all 
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the states are localized at certain coordinates X^, these 
transitions correspond to "jumps" of a particle in the 
X coordinate space from the position Xi^ to positions 

Xk- Therefore the operator L^'^ corresponds to an in- 
stant jump with the jump lengths distribution due to 
the exponential law in accordance with the overlapping 
integrals Ai{Xk). Another operator Aqq^ changes only 

the phase of the complex amplitude . Time duration 
of the action of this operator is tfj/n, which is different 
for different powers n. Now we introduce a probabil- 
ity distribution function (pdf) V{X,t) = \Ck{t)\^ to be 
a particle at position X at time t. Since the dynamics 
of the pdf |CA;(i)P is determined by the same Liouville 
operator as in Eq. namely 



dt\Ckit)\^ ^ L\Ck{t)\'' , fc=l,2,. 



(21) 



we obtain that the pdf corresponds to the exponen- 
tially localized solution [H V{X, t) ~ A\{X) sin^ (Mo* 

which is relevant to Eq. (fT9|) for t < tp. 

For t ^ tp Eq. (PT|) we present in the integral form 



riX,t)= / L\Ckit')\'dt' 

JQ 



(22) 



Since, in the new terminology, summation over indexes 
k corresponds to integration in space, the r.h.s. of Eq. 
can be rewritten in the form of the integral operator 



dX'V{X,t]X't')V{X',t' 



t pt 
L\Ck{t')fdt' / dt' 

^0 



In this case all combinations of the overlapping integrals 
A(h) with corresponding differentiating over Uk play the 
role of the kernel or transition probability 'P{X,t; X't') 
of this transformation. Therefore one has to consider a 
variety of combinations of the operators 



r. 9 1 


pi 
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92 




P2 




L 90/ J 












L d(j>ki 







where pi + q2 + ■ ■ ■ + qi + P2 + ■ ■ ■ ~ P + Q — n. This 
corresponds to different realizations of instant jumps and 
waiting times between any two successive jumps. Note 



that 



Ao-: 



pi 



corresponds to waiting time of duration 

Pit/ P. To proceed, we follow ideas of so-called continu- 
ous time random walk [2^[3^. Therefore, the transition 
probability consists of the pdf of jump lengths f{X — X') 
and the pdf of waiting times w{t — t'). For simplicity, we 
suppose that V{X, t; X't') = f {X - X')w{t - t'). 

From the exponential decay of the overlapping inte- 
grals on the large scale one obtains that all jump lengths 
have finite expectation values and variances. Note, that 
" " (0)1 r ' '■ (0)1 "^^ 

L\' Li — L\ ' corresponds to a "jump" which is 

a composition of random walks (e.g., the simplest real- 
ization is presented in [131) • For large Qi the displace- 
A 

tral limit theorem. Here A; = Xi 



ment '-^i Gaussian distribution due to the cen 



X'l are transition 



lengths due to operator either L^"-* or Li — Lf K There- 
fore, we can believe that these lengths also are approx- 
imately distributed by the Gaussian law f{X — X') = 
exp(-A2/2CT^)/V27rCT2 and = (A^). It is worth stress- 
ing that the overlapping integrals do not specify the pdf 
of waiting times. Therefore, ^(t) can be defined from the 
average value of the waiting times (r) = Tw{T)dT. 
This value also can be calculated from the following ar- 
guments. For asymptotically large i, waiting times are 
T = pt/P, where p G [1,-P] and P G Therefore, 
one obtains for the average waiting time 

Tl p 

(r)= Km ly±ypi = t. 

^ ' n->oo n P ^ P 

p=i p=i 

This value diverges with t — > oo, and it means that there 
are infinitely many realizations of waiting times of the 
order of t [31|. To fulfill this condition, the waiting times 
are distributed by power law w{t) — atp/r^'^"' where 
< a < 1, such that (see [sH ) 



at^ 

T-l+a 



dr — I and 



rdr 



oo . 



(23) 



It is reasonable to suppose that random jumps and wait- 
ing times are independent and identically distributed pro- 
cesses. Therefore this random qualitative description of 
Markov operator L in Eq. pT|) and Eq. respectively, 
corresponds to the continuous time random walk (see e.g. 
Refs. [H, m, [l^l) which is described by the fractional 
Fokker-Planck equation 



dtV{X, t) - Do.dl-"dj,V{X, t) = , 



(24) 



where Da = c^/t^ is a generalized diffusion coefficient 
and is a designation of the Riemann-Liouville frac- 
tional derivative 



dtfit) 



1 



r(-t.)yo (t-T) 



f{r)dT 



In this case only nonzero initial conditions can be taken 
into consideration. Without restriction of the generality 
one can consider 'P{X,t = 0) = S{X ~ Xig). Equation 
((24| describes subdiffusion [H, [33, 34] since a < 1. From 
Eq. ([M)) one obtains for the second moment 



(X'it)) 



X^ViX,t)dX = Dat" 



(25) 



The transport exponent a cannot be specified here from 
the developed probabilistic arguments. In the recent nu- 
merical studies of the discrete NLSE [ll,[li| the exponent 
a was found in the range 0.3 0.4. 

In conclusion, the dynamics of the initially localized 
wave packet il){x^t) was studied in the framework of 
the NLSE in the presence of a random potential. It is 
shown that the dynamics may possibly be described in 
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the framework the Liouville operator. The interplay be- 
tween disorder and nonhnearity /J leads to the compli- 
cated dynamics of the initially localized state ip{x,t — 
0) = 'i'ui„{x). So, the influence of the nonlinearity on the 
initial time dynamics is weak, and a perturbation theory 
in (3 can be developed. An analytical expression for a 
wave function of the initial time dynamics is found by 
the perturbation approach. As follows from a perturba- 
tive solution, at the initial times t < 1/ (3 the nonlinearity 
affects mainly the phase of the wave function, while the 
shape of the wave packet corresponds to the exponential 
localization due to the overlapping integrals described by 
Eq. dUl). 

At asymptotically large times t ':$> 1/(3 the nonlinear 
effects become important. To evaluate the influence of 



the nonlinearity on the rate of spreading of the initial 
wave packet, one can consider the large times asymptotic 
dynamics of the tails of the packet. In this case, the 
dynamics can be described qualitatively in the framework 
of a phenomenological probabilistic approach by means 
of a probability distribution function V{X,t) = |Cfe(t)p. 
The last may be governed by the fractional Fokker-Planck 
equation which describes the asymptotic behavior of 
the tails of the wave packet; its solution corresponds to 
subdiffusive spread of the initially localized wave packet. 
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